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Abstract. New methods of solutions of the DGLAP equation and their implementation through 
NNLO in QCD are briefly reviewed. We organize the perturbative expansion that describes in x- 
space the evolved parton distributions in terms of scale invariant functions, which are determined 
recursively, and logarithms of the ratio of the running couplings at the initial and final evolution 
scales. Resummed solutions are constructed within the same approach and involve logarithms of 
more complex functions, which are given in the non-singlet case. Differences in the evolution 
schemes are shown to be numerically sizeable and intrinsic to perturbation theory. We illustrate 
these points in the case of Drell-Yan lepton pair production near the Z resonance, analysis that can 
be extended to searches of extra Z'. We show that the reduction of the NNLO cross section compared 
to the NLO prediction may be attributed to the NNLO evolution. Q 
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INTRODUCTION 

QCD will enter its precision era with the advent of the LHC. For this reason, from the 
theory side, precise determinations of some key observables at LHC energies will be 
essential both for the discovery of new physics and for QCD partonometry. One of the 
first searches that we expect to be performed at the LHC will be the study of the rapidity 
and the invariant mass distributions in Drell-Yan on the Z resonance and at larger invariant 
mass values, which will allow to calibrate the parton distributions functions (pdf's) at the 
new (higher) factorization scales available at the new collider. In return, these first studies 
of the background will be crucial in order to quantify the luminosity of the machine and 
for the investigation of possible extensions of the Standard Model, from the "simplest" 
ones concerning the Higgs or extra neutral gauge interactions, to more complex ones 
involving supersymmetry or popular scenarios such as the those concerning theories with 
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extra dimensions. Given the high accuracy of the current computations, which in some case 
have reached the few-percent level of precision, and the rather small variations observed 
in moving from next-to-leading order (NLO) to next-to-next-to-leading order (NNLO) in 
some key observables, it is rather natural to ask whether the sistematic errors due to the 
implementation of the evolution of the pdf 's are of the same order of the new perturbative 
contributions coming from the hard scatterings. We answer these points by presenting 
a new method of solutions of the DGLAP equation that we have developed and tested 
in several previous works that we call the method of the logarithmic expansions. The x- 
space solutions are classified in terms of the different types of accuracies selected in the 
construction of the evolved pdf 's, solutions that bring in a systematic uncertainty in the 
prediction of the hadronic observables. The theory behind this new approach can be found 
in [1,0]. 



THE NNLO EVOLUTION 

Studies of hard scattering cross sections such as the invariant mass distribution of the 
lepton pair in Drell-Yan do /dM [|3|], and the rapidity distribution do/dMdY f?] are 
available and can be used as tests for the NNLO evolution. Some of these computation 
have been performed before that the analytical NNLO DGLAP were computed (see [5]). 

We recall that the perturbative expansion of the DGLAP splitting functions and of the 
/3 -function take the generic form (to the m-th order) 



with /3yt being the corresponding coefficients of the /3 function which have been summa- 
rized in llj]. Leading order (LO), NLO and NNLO expansions correspond to the cases 
m = 0, 1, 2 respectively. 

We solve the DGLAP equation by formulating a general ansatz in (Bjorken) .x:-space 
which holds to any perturbative order and which allows us to have a control on the accuracy 
in terms of powers of the strong coupling constant. It is given by ^ 
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The expansion is valid both in the non-singlet and in the singlet sectors, and the scale- 
invariant coefficients A^{x) are determined by solving a chain of recursion relations ob- 



tained by substituting the ansatz into the DGLAP equation. We have implemented the 
method in a numerical code Candia, which will be made available soon. 



The logarithmic resummation of the DGLAP solutions 

When we move from logarithmic expansions of the form (O to exact solutions, we are 
actually performing a resummation of theperturbative expansion of the pdf's. This point 
has been discussed in our original works [lli|2|], to which we refer for more details. Here 
we just quote the result. 

The general solution can be written in terms of A'j,5'^,C^, coefficients that will be 
calculated by a chain of recursion relations PJ] giving 

f(^,<f) = fi^^") it'^^A (f^A f(-.Qi) 
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where ® is the standard convolution product of the DGLAP and where 

DlXx) =A'„{x)®B[_„{x)®C',_,{x)®f{x,Ql). (4) 
The explicit expression of the functions appearing in the ansatz are 

^ = log^, (5) 
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and 
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S£ now corresponds to the LO solution, which is the same as in (O while M and ^ 
correspond to new functions of the running couplings at the two evolution scales, initial 
(2o) and final {Q). These two functions perform a resummation of the simpler logarithms 
included in At NLO the implicit resummation performed by the exact non-singlet 
solution can be understood from its explicit form in x-space 
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which generates the ordinary contributions log f ^ ) after an expansion 
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LHC PHENOMENOLOGY 



Both in the evolution performed with Candia and in the MRST and Alekhin evolution 
[|6i|7|] we use the same inputs given by these authors, choosing the initial scale /ig = 1.25 
GeV^, and the same treatment of the heavy flavors, a rather nontrivial point from the 
numerical side. On the Z resonance we get for the NNLO ^-factors the values 

K{M-£) = (a^wL0®*M^5?)/(^^L0®*M^?r) =0.97 
K{Mz) = (^^^LO®*SfNDiA)/(^^^o®*CA%A) =0-95 

K{Mz) = (O^NLO^^TkhJ/i^NLO^^Tklun) =0.9S (12) 

where 4> denote the parton luminosities, which correspond to a reduction by 2.7% of the 
NNLO cross section compared to the NLO result (MRST evolution) and larger for the 
Candia evolution (4.4%), while for Alekhin is 1.5%. From the analysis of the errors 
on the pdf's to NNLO, for instance for the Alekhin's set, the differences among these 
determinations render the results compatible, being the variations on the ^-factors of the 
order of 4% (see the discussion in [20). 

One can estimate the role played by the evolution, respect to the NNLO hard scatterings 
by separating the two contributions in the variation of the hadronic cross section Ao 
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FIGURE 1. Cross Sections in the region of the peak of the Z boson at LO, NLO, and NNLO obtained 
using the luminosities evolved respectively by MRST 
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FIGURE 2. Cross Sections in the region of the Z with a zoom in the peak region. 



between the NNLO and the NLO determinations, into contribution from the change of 
the hard scatterings (Ad) and from the luminosities (A0) 

A(7 = Aa®(^ + a®A(^. (13) 

A numerical analysis shows that the change in the luminosities play a dominant role in the 
overall variation and_are induced by the new terms in the NNLO evolution. Some of our 
results, taken from [2] are shown in Figs. [T] and [21 The NNLO reduction for the invariant 
mass distributions is quite evident respect to the NLO result. 



CONCLUSIONS 



In the cross sections that we have studied the NNLO QCD corrections are small (at 
the few percent level), but the reduction of the dependences from all the scales is 
significant. At the same time the issue of how we solve the evolutions equations is 
not secondary, given the small variations involved and the claimed accuracy. Given 
this issue and the problem with this systematic source of errors, our confidence on 
the parameterizations of the pdf's should never be too high, and it is reasonable to 
keep a critical edge. We recall that with 10~^ fb of integrated luminosity the sta- 
tistical error expected on the Z peak is around 0.05% at the LHC, much smaller 
than the systematic error intrinsic to an NNLO prediction due to the evolution. 
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